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Abstract. The spectral evolution of powerful double radio galaxies (FR IPs) is thought to be determined by the 
acceleration of electrons at the termination shock of the jet, their transport through the bright head region into 
the lobes and the production of the radio emission by synchrotron radiation in the lobes. Models presented to 
date incorporate some of these processes in prescribing the electron distribution which enters the lobes. We have 
extended these models to include a description of electron acceleration at the relativistic termination shock and 
a selection of transport models for the head region. These are coupled to the evolution of the electron spectrum 
in the lobes under the influence of losses due to adiabatic expansion, by inverse Compton scattering on the 
cosmic background radiation and by synchrotron radiation. The evolutionary tracks predicted by this model are 
compared to observation using the power/source-size (P-D) diagram. We find that the simplest scenario, in which 
accelerated particles suffer adiabatic losses in the head region which become more severe as the source expands 
produces P-D-tracks which conflict with observation, because the power is predicted to decline too steeply with 
increasing size. Agreement with observation can be found by assuming that adiabatic losses are compensated 
during transport between the termination shock and the lobe by a re-acceleration process distributed throughout 
the head region. 
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1. Introduction 
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Powerful double radio galaxies or 'classical doubles' (CDs) 
owe their name to the extended (hundreds of kpc) lobes of 
radio emission they exhibit on opposite sides of the par- 
ent galaxy. Fanarofffe Riley (1974| ) classified these sources 
as type II objects; they have luminosities Pi78MHz > 
5 x 10 25 WHz -1 and are edge-brightened, with bright 
outer hotspots. It is universally agreed that the radio con- 
tinuum of the CDs is synchrotron radiation from relativis- 
tic electrons and perhaps positrons. The standard scenario 
is that of a jet propagating from the galaxy to the outer 
parts of the lobes, passing through a shock front at the 



of the galaxy. Subsequently, Kaiser fe Alexander (1997 ) 
showed that the cocoons also expand self-similarly. Based 



on this result, Kaiser et al (1997) developed an analyti- 
cal model for the spectral evolution of FR II sources as a 
function of redshift, jet power and the scaling of the exter- 
nal density profile. They assumed that electrons with an 
initial power-law distribution in energy are continuously 
injected into the plasma immediately downstream of the 
shock front that terminates the jet. A key point of their 
treatment is that the pressure in this region equals that 
which drives the bow shock into the external medium on 
the axis of the source. This results from the fact that the 



hotspo ts and subsequently filling a 'cocoon' around the 



jet is confined by the pressure of the cocoon (Begelman 



jet wit h radiating particles see, for example, Begelman 



fe Cioffi (1989); Peacock (1999) 



Cioffi, 1989), so that the working surface over which its 



Faue (1991) considered radio sources as expanding 
bubbles that are fed by supersonic jets and drive a 
bow shock into a radially stratified external medium. He 
showed that the jet length and the bow shock grow in 
a self-similar way for external atmospheres in which the 
density drops off more slowly than 1/r 2 from the centre 
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thrust is distributed expands self-similarly, along with the 
bow shock and cocoon. Particles undergo adiabatic losses 
on moving from the terminal shock front into the main 
part of the cocoon, whose pressure is lower than that at 
the working surface by a constant factor (taken as rs 16 
by Kaiser et al (1997 )). They then undergo synchrotron 
and inverse Compton losses, as well as adiabatic lo sses as 
the lobe of the source expands. Using this model, Kaiser 
et al (1997| ) computed evolutionary tracks in the power- 
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University of linear size (P-D) plane. These display decreasing power as 
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the size increases, in agreement with observations, which 
indicate a deficit of large luminous sources. 



Bow shock 



However, Blundcll et al (1999) noted that there is no 
evidence that the size of the working surface is propor- 
tional to the size of the source. Instead, the hot-spots of 
all known FR II sources appear to be a few kiloparsecs in 
diameter. This has a profound influence on the predicted 
evolutionary tracks, because it implies that the adiabatic 
losses suffered upon moving into the cocoon by particles 
accelerated at the termination shock are not constant, but 



increase strongly with the age of the source, plundcll et al 
(1999)]] modelled the spectral evolution of CDs including 
this effect by prescribing the spectrum of electrons en- 
tering the cocoon to be a broken power-law distribution. 
They, too, found evolutionary tracks that agree with the 
observed lack of large, luminous radio sources, as well as 
with several other properties of the samples they investi- 
gated. 

In this paper we present computations of evolution- 



ary tracks based on the picture adopted by Blundcll et al 
(1999|]] However, instead of prescribing the electron distr 
bution that enters the cocoon, we assume that electrons 
are accelerated by the first-order Fermi process at the ter- 
mination shock and then propagate through the hot spot 
region into the cocoon or 'lobe' according to one of two 
models, which we designate case A and case B. In case A, 
the full adiabatic energy loss corresponding to the age 
dependent pressure difference between the hot spot and 
cocoon is applied. This results in P-D tracks which are 
in conflict with the observations. One possible way out 
of this problem is to assume that a re-acceleration pro- 
cess occurs after the initial encounter with the termina- 
tion shock. Such a process is indicated independently by 
investigations of the spectra of indiv idual hot-spots and 



of o ptical s y nchrotron emitting jets ( Meiscnheimer et al 
1996|, Il997fc |Perley et al, 1997t |Wagner fc Krawczynski, 



2000]^ This motivates our case B model, in which we as- 
sume that the adiabatic losses between the termination 
shock and the lobe are compensated by a re-acceleration 
process during propagation through the head region. In 
both cases electrons are carried through the high-loss re- 
gion by fluid elements whose residence time is distributed 
according to a specified transport equation. During this 
time they suffer synchrotron and inverse Compton losses. 
After entering the cocoon, we follow the particle distri- 
bution as it cools and radiates and compute evolutionary 
tracks as well as snapshots of the spectrum at different 
ages. 

The paper is organised as follows: in Sect. 2 we sum- 
marise the hydrodynamic picture which we adopt for the 
source evolution. In Sect. 3 our treatment of Fermi accel- 
eration is described, the transport model is specified and 
the kinetic equation obeyed by particles in the cocoon is 
formulated. Section 4 is devoted to a series of tests and 
special cases to illustrate the way in which the transport 
and acceleration models influence the particle distribution 
in the hot spot and the lobe. Our main results on spec- 
tral evolution are presented in Sect. 5, where we conclude 




Contact discontinuity 



Shocked intergalactic medium 
Hot spot 
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Fig. 1. A cartoon of the interaction of the intergalactic 
medium and a jet of an FR II radio galaxy, showing the 
principal components of the spectral evolution model. A 
jet of constant Lorentz factor and power is decelerated in a 
hot spot (whose size is independent of source age) convert- 
ing part of its energy into relativistic electrons and part 
into magnetic field. The electrons propagate according to a 
prescribed transport model [parameterised by a and r — 
see Eq. (|l8|)] from the hotspot into the turbulent head 
region. The hot spot moves around inside the head re- 
gion; averaged over time, the jet thrust is balanced by the 
pressure in the head region, which expands self-similarly 
along with the bow shock and lobe. Synchrotron losses 
are important in the (constant) magnetic field of the hot 
spot, but adiabatic losses suffered in moving into the head 
region and later into the lobe may be compensated by a 
re-acceleration process acting in either the hot spot or 
the head region. Once in the lobe, time-dependent syn- 
chrotron, inverse Compton and adiabatic losses determine 
the spatially integrated electron spectrum. 



that the observed data exclude case A i.e., they can be ex- 
plained only by a model such as case B, that includes re- 
acceleration of electrons after their encounter with the ter- 
mination shock. This conclusion and our results on source 
spectra at different ages are discussed and compared with 
previous work in Sect. 6. 



2. Hydrodynamic Model 

The standard model for the CD radio sources comprises 
a central object (AGN) from which emanate two jets in 
opposite directions, embedded in a cocoon or lobe of dif- 
fuse emission. The jets themselves terminate in regions 
of intense emission, called the hotspots or working sur- 
faces, after which the shocked jet material flows into the 
cocoon. This material drives a bow shock into the external 
medium, from which it is separated by a contact disconti- 
nuity. A cartoon of the model is shown in Fig. |l|. 

Following iBlundell et al (199S| ), we assume that the 
radial profile of the external (intergalactic) medium does 
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not vary with redshift and scales with distance r from the 
centre of the host galaxy as 



Pc{r/r c )~ 



(1) 



flGarrington fc Conway, 1991 ; Mulchaey fc Zabludoff. 



1998|) . |Typical values for the above parameters are: p c = 



1.67 x !CT 23 kg m 



-3 



= lOkpc and (3 = 1.5. 



During the lifetime of the source, we assume that the 
power Qq carried by the jet and its Lorentz factor rj et re- 
main constant. The jet terminates at a shock front which 
we identify as the primary hotspot embedded in a high- 
pressure region at the end of the lobe which, following the 
terminology of Blundcll ct al (1999 ), we call the head re- 
gion. The exact position of the termination shock within 
this region is assumed to vary in the manner proposed 
by Schcucr (1982) in the 'dentist's drill model'. Averaged 
over time, the thrust of the jet is distributed over the head 
region, which itself expands with the rest of the source. 
However, the size of the primary hotspot is independent 
of source age. Although its ram pressure substantially ex- 
ceeds the average pressure in the head, simulations indi- 
cate that it does not advance much more rapidly (Cox et 
al, 199[[), because it quickly encounters material of higher 
than average density. Consequently, the pressure ph im- 
mediately downstream of the termination shock can be 
estimated as the jet thrust, Qo/c, divided by the area, A, 
of the primary hotspot itself. Assuming the jet is relativis- 
tic and the magnetic field is amplified to reach equiparti- 
tion with the thermal energy density, we may express the 
(constant) magnetic energy density in the hotspot, ith, in 
terms of the jet power as follows: 



(2) 



3Qo 
2 r " 2cA 

' Falle (19 9?) and |Komissarov k Falle (1998j ), found so- 
lutions in which the size of a radio source increases in a 
self-similar way with time, with the total length scaling 
as: 



D(t) = 2ci 



t 3 Qo 



1/(5-/3) 



(3) 



From the measured expansion speeds of powerful sources 
one concludes that c\ « 1.8 ( Blundcll ct al, 1999 ). 

In contrast to the constant volume and magnetic field 
of the hotspots, the lobes of the source expand self- 
similarly. Applying the jump conditions at the external 
bow shock and using Eq. (||) , gives for the pressure of the 
downstream plasma immediately behind the bow shock: 



Pd{t) = 



18 



,(2-/3) 



r x + 1 (5 - /?) 2 



(per?) 



0\ 3/(5-/9) 



y(2-j3)/(5-l3) t -(4:+/3)/(5-0) 



(4) 



( Kaiser fc Alexander, 1997 ) where T x is the ratio of specific 
heats of the external gas, which we take to be 5/3. 

The pressure in the lobes of the source, piobe, is some- 
what lower than this value, depending on the axial ratio 



of the source. In this way the pressure gradient is provided 
to drive a high speed back-flow from the head towards the 
core of the galaxy (Liu et al, 1992). Assuming, in addition, 
that the magnetic field is tangled and that equipartition 
is also attained in the lobes, leads to the conclusion that 
the magnetic energy density in the lobe, uiobcj scales as: 



U\obe(t) = Uh 



(5) 



where a — (4 + /3)/(5 — f3) and to is the time at which 
the size of the head region was comparable to that of the 
primary hot spot: 

V« /„ _ B\ 3/(4+/3) 



tn = 



2-/3 



r x + i (5 - p) 



:CA 



(6) 



3. Particle acceleration and transport 



We assume that the first-order Fermi process operates at 
the termination shock in the jet and accelerates there elec- 
trons (and possibly positrons) into a distribution which 
can be described by a power-law in energy. As they move 
away from the shock these electrons radiate and form the 
primary hot spot within the head region. We can consider 
this process as an injection of electrons at a rate Q s (j) 
into the plasma behind the termination shock: 



s(7) 



g 7 P for 7min < 7 < In 



otherwise. 



(7) 



If the jet Lorentz fact or Fj ct is of the order of 10, we can 



choos e 2 < p < 2.3 ( [Kirk et al, 2000| ; [Achterberg et al 
|200 1| ). The upper and lower limits on the Lorentz factor 
7max and 7 m i n are not available from theory. It is plausi- 
ble to assume 7 m i n « Tjot , but 7 max should be determined 
by a local balance of the acceleration and cooling rates. 
Fortunately, our models are not sensitive to this parame- 
ter, since synchrotron cooling within the head region will 
in any case determine the effective maximum energy of 
those particles which enter the lobe. The fraction of the 
jet power, rjQ , (with rj < 1) which is transfered into accel- 
erated particles at the termination shock is J ^m c c 2 Q s &y. 
By means of Eq. (Q) this gives for qq: 

vQo 



Qo = 



m e c 



r(p-2)( 7) 



2-p 
min 



/ max 7 



(8) 



After injection into the primary hot spot, electrons 
are transported by the turbulent motions of the plasma 
through the head, and eventually reach the lobes. During 
this time, which we assume is short compared to the life- 
time of the source, they are subject to synchrotron losses 
in the relatively strong magnetic field generated behind 
the shock as well as inverse Compton losses by scattering 
off the cosmic background radiation. The effect which this 
has on the power-law distribution injected at the termina- 
tion shock depends on the distribution of the 'escape time' 
i.e., on the probability an injected particle has of leaving 
the region of strong field after a certain time. 
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This distribution can be modelled in several ways. The 
simplest approach would be to assume each particle es- 
capes after a fixed time, corresponding, for example, to 
the time taken for a fluid element to cross the head region 
at a given speed. Upon reaching the lobes, the particles 
would then have the distribution computed by Kardashev 
(1962JJJ which vanishes above a cut-off energy at which the 
escape time equals the cooling time. However, the turbu- 
lent appearance of the hot spot and head regions indi- 
cates that not all elements of shocked plasma will take 
the same time to escape into the lobes. At the next level 
of sophistication, this can be modelled in a diffusion pic- 
ture — a plasma blob, and along with it accelerated par- 
ticles, performs a random walk through the head region. 
A convenient formalism, which has been widely used to 
describ e transport in turbulent flui ds and other random 
media ( [Bouchaud fc Georges, 1990 ) is giv en by the the- 
ory of continuous time random walks (see Ragot & Kirk 
(1997|]]and refer ences therein). Diffusion in space results 
if a sequence of jumps of fixed length a in a random direc- 
tion is performed, where each jump takes the same time 
td • The advantage of this formalism is that by changing 
the waiting time distribution between jumps it is easily 
generalised to include non-diffusive or 'anomalous' trans- 
port characterised by a single transport parameter a. The 
mean square distance (Ar 2 ) of a particle from its starting 
point is linearly proportional to the elapsed time only in 
the special case of diffusive transport a = 1; in general one 
has (Ar 2 ) oc t a with < a < 2. Sub-diffusive transport 
(a < 1) occurs when the motion of the fluid elements is 
confined in some way — for example to forwards and back- 
wards motion in a channel, which itself diffuses sideways. 
Supra-diffusive transport, on the other hand, corresponds 
to a more 'free-streaming'-like motion, less impeded by the 
random eddies. The spatial propagator for the generalised 
transport process in spherical geometry is 



P(r,t) 



3 ^OToy ^2(1-^3/12(2^)] ^2 




— a 



3r 2 /aT D \ a 
2^ \~TJ 



ct 3 tt 3 / 2 
-,i/(2- Q n 



(9) 



where td is a characteristic timescale for the jumps. 

During propagation, the rate at which an electron of 
energy j(t)m e c 2 loses energy is described by the equation: 



d", 
df 



37 = -57 



fJ 



4 cr T , 
3 m c c 



(10) 



where ctt, fn c and c are the Thomson cross section, the 
electron rest mass and the speed of light, respectively. The 
energy density u c of the cosmic background radiation field 
is given in terms of an effective magnetic field B c by 



2^0 ' 



B c = 0.318(1 



z] 2 nT 



(11) 



and is taken to be constant during the lifetime of a radio 
source. Equation ( [Toj ) may be integrated to give r )~ 1 (t) = 



7~ 1 (£ / ) + (t — t')g. These loss processes may be included 
on the right-hand side of Eq. (^J) , leading to a composite 
propagator P(x,t, 7,7') that expresses the probability of 
finding a particle with Lorentz factor 7 at x and t, given 
that it started at the origin (x = 0, t = 0) with Lorentz 
factor 7': 



P(x,t, 7 , 7 ') = P(x,t)5 



7 . 



7 



1 + tgi 



(12) 



The distribution Nr(^) in Lorentz factor of the total 
number of particles contained in a sphere of radius R, 
at the centre of which particles are injected according to 
Eq. (M) starting at time t stai t is 



pR />00 r-t 

N R (j) = 4tt / drr 2 / d 7 ' / dt' 
P(r,t-t', 7 , 7 ')Q s (7') ■ 



(13) 



The two timescales that determine Nr are the transport 
timescale (a generalisation of the diffusion time) 



t diS = (2 - a )V-«)/° 



3(2 - a)R 2 
2^ 



l/a 



CVT D 



(14) 



and the cooling time, defined for a particle of Lorentz fac- 
tor 7 m i n as icooi = l/(<77min)- If these are both short com- 
pared to the age of the source, it is reasonable to assume 
the transport process into the lobes reaches a steady state, 
found by setting t s tart — > — 00: 

'0 for 7 > 7, nax 

N R(l) = { for 7min < 7 < 7max (15) 



for 7 < 7 min 



where s = (7' — 7)/ (577')- The kinetic equation satisfied 
by Nr reads 



'''- V " ^ ( 91 2 N R ) = 0,( 7 ) - Q h (7) 



Of 



97 



(16) 



where Qh(7) is the rate at which particles leave the sphere 
of radius R (and so enter the lobe). In the steady state, 
it is straightforward to solve Eq. (|l^) for Qh- The result 
is conveniently expressed as a modulation function M de- 
scribing the effects of cooling on the incident distribution 
Q s , and defined as 



Qh(7) = 9o7 p M{JiP<i,P,t) 



(17) 



where we have also defined parameters representing the 
Lorentz factor: 7 = 7/7 m i n , the reciprocal of the dynamic 
range of the incident distribution: pd = 7 ro in/7max and the 
dimensionless ratio of the transport time to the cooling 
time of a particle at 7 m i n : r = fdiff Acooi- From Eq. (16) 
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wc find 



A^(7,Pd,P,T) 



for 7 > l/p d 



4= f^dxa; 1 ^ 
Jx i 



1 



= < 



for 1 < 7 < l/pd 



p-2 



(18) 



where 



X2 = 



7r 



(1 - Pdl) 



for 7 < 1 



a/(2-a) 



r 7 



P-2 



7T 



Note that transport from the primary hot spot into 
the lobes results from the diffusion of plasma elements in 
a turbulent environment. The mean free paths involved 
are much larger than those estimated for individual ener- 
getic electrons, which are usually taken to be of the order 
of the gyro radius. The turbulent mean free path is also 
independent of the particle energy. However, this expres- 
sion still does not account either for adiabatic losses or for 
any acceleration process that might occur between leaving 
the termination shock and entering the lobe. 

The effect of adiabatic losses is to reduce the energy of 
each particle by a factor given by the ratio of the pressures 



before and after the expansion [see, for example, Scheuer 
& Will iams (1968| )]. In the case of a transition from the 
pressure in the primary hot spot to the pressure in the 
lobe we have 7 — > (ui b c /uh) 1 ^ 4 7- As a result, the slope of 
the power-law part of the distribution remains unchanged, 
but the upper and lower limits shift downwards in energy. 
Of course, the situation is more complicated if, simulta- 
neously, radiation losses are important. Our treatment of 
losses in the lobes, where an explicit time dependence of 
the energy density is available, takes this explicitly into 
account. Qualitatively, however, the form of the spectrum 
should not depend on the sequence in which these pro- 
cesses operate. In treating propagation through the head 
region we impose the adiabatic losses after the radiation 
losses, bearing in mind that the values derived for quanti- 
ties such as the diffusion time idiff may depend somewhat 
on this choice. Thus, assuming constant injection at the 
termination shock and requiring particle number conser- 
vation, we find for the function describing the injection of 
electrons into the lobes: 



Hob. 



,(7,t) = kQ h (kj,t) 



(20) 



where k(t) = [uh/uiobc(i)] 1/4 - 

Re-acceleration by multiple encounters with shock 



fronts has been considered in many papers, notably Spruit 
1988|)]|Achtcrbcrg (1990|); |Anastasiadis fc Vlahos (1993|) 



Schneider (199$; Pope fe Melrose (199J ); |Marcowith fc 
Kirk (t999| ) ; |Gieseler fe Jones (2000|) . In our model, we 



distinguish between the termination shock that forms the 
primary hot spot and a sequence of subsequent, weaker 
shocks through which the plasma passes before emerging 
into the relatively inactive lobe region. The theory of dif- 
fusive shock acceleration predicts that a single encounter 
with a weak shock will not change the power-law index 
of an incoming distribution produced at a strong shock 
- the effect is merely to increase the amplitude of the 
distribution, and to produce a steep power-law tail above 
the upper cut-off energy of the incoming particles. Passage 
through many such shocks, however, tends also to harden 
a preexisting power-law if the ratio of the escape time to 
the time between shock encounters is large. Shock drift 
acceleration, on the other hand, produces an effect which 



is similar to an adiabatic compressio n ( Begelman fc Kirk 



fiqj 1990 ; Anastasiadis fc Vlahos, 1993 ), in that the power- 



law index is unchanged, but the distribution is shifted and 
slightly modified around the cut-off energies. The effect 
of a highly turbulent environment comprising of a 
mixture of different types of accele rators has been 
studied by Manolakou et al (1999 ) who conclude 
that the level of turbulent activity has at most a 
minor impact on a pre-existing power-law index, 
and contributes only to the efficiency of the pro- 
cess. Given the variety of possible effects, it is clearly 
a difficult observational challenge to distinguish between 
the different mechanisms. In this paper, we do not adopt a 
specific acceleration model, but simply assume that when 
re-acceleration is important a distribution emerges into 
the lobe which can be described as a power-law above a 
lower cut-off energy, modified by synchrotron and inverse 
Compton losses at high energy. This corresponds precisely 
to a modification of the effectiveness of the adiabatic losses 
described by the parameter k(t). If we now interpret the 
quantity rj as that fraction of the jet power which goes into 
energetic particles which enter the lobe, then the effects 
of re-acceleration correspond simply to setting k{t) = 1 in 
Eq. ©. 

3.1. The evolution of the electron distribution in the 
lobes 

The spatially integrated distribution N(~f,t) of particles 
in the lobes, is governed by a kinetic equation similar to 
Eq. ©: 

ON d 



dt dj 



UN) =Qlobe(7,*) 



(21) 



However, in this case the loss term is more complicated, 
involving adiabatic losses as well as synchrotron and 
Compton losses, the first two of which have an explicit 
time-dependence: 



7 



.7 



6ic7 2 



11 

5 r 



(22) 



The three terms on the right hand-side of Eq. (|22j) rep- 
resent (i) adiabatic losses, with b = a/4 according to 
Eq. (|^), (ii) losses by inverse Compton scattering, with 
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&i c = (4/3)(crT/mec)w c a constant, and (iii) synchrotron 
losses in the adiabatically decreasing magnetic field of the 
lobe [see eq. (§)]: b s = (4/3)(<r T /m c c)u h t%. 

Following the method of characteristics, Eq. (|2l]) may 
be rewritten in terms of the derivative of N along a fam- 
ily of curves x(l->t) describing the evolution of individual 
particles and given by integrating Eq. (E3): 



1 

7> 



b\ c t 



bd 1 



(23) 



1-6 l-b-a 
Then N(x,t) satisfies the ordinary differential equation 



dN 
~dt 



+ 2j(b ic + b s t- a ) 



N = Qlobe(7j*)> 



(24) 



in which 7 is considered a function of \ and t according 
to Eq. ( |23|) 

Eq. (p4|), with Qiobe given by Eq. (|20|), may be in- 
tegrated numerically, the upper and lower limits of inte- 
gration being the observation time, t and the earliest time 
ti (-7, t) at which the source injected particles into the lobes 
that subsequently cooled to the chosen Lorentz factor 7 
at time t. The initial condition is N(x,ti) = 0. 

The injection time t\ is a function of the Lorentz factor, 
7 and the observing time, t. It is determined not only by 
the cooling rate, but also, at least potentially, by the time 
^start at which the source started its activity. This limit is 
important if particles injected into the lobe at time i s t a rt 
with 71 < 7max/fc(^start) can cool to 7 at time t. From 
1), this situation arises when 



Eq 
1 





- = ts 



Ms 



+ t s 



h xii,t) 



< 



k(t s tait) 



(25) 



7max 

in which case we set t\ = £ s tart- If, on the other hand, 
7i > 7max/fc(istart), then ij, is found from the relation: 

k(t) . ( h r Kt7 a 



On 



1 



1 



(26) 



In summary, the injection time, t\ lies in the range 
tstart < U < t, if 7; > 7 m a X /M i start) and equals t start if 

7i ^ Tmax 

Finally, the specific power P v (t) (power per unit fre- 
quency) is found by integrating the product of AT (7, t) and 
the emissivity of a single electron. Adopting the delta- 
function approximation we have: 



1 f c 



ao7 2 SLW 2 <5(^-ai J B L (t)7 2 )Ar(7,i)d 7 (27) 



where B^,(t) is the magnetic field in the lobe at time t, 
a Q = 1.6 x 1(T 14 WT- 2 and ai = 1.3 x lO^HzT" 1 . 

4. Tests and special cases 

4.1. Cooling in the head region 

The effect of cooling in the relatively strong magnetic field 
of the head region is described by the modulation function 
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Fig. 2. The modulation function M(j, Pd,P, t) defined in 
Eq. ( jig ) as a function of 7, for pd = 10~ 4 , p = 2.23 and 
for a = 1 (diffusion), a — 0.5 (sub-diffusion) and a = 1.5 
(supra-diffusion) . In each case the solid, dotted, dashed 
and long-dashed curves correspond to r = 10~ 4 , 1CP 2 , 1 
and 10 2 , respectively. 



M of Eq. @. This function is plotted in Fig. |for various 
values of the ratio of diffusion to cooling time r. 

The character of the three different transport regimes 
is clearly seen from this figure. In the case of supra- 
diffusion, a > 1, the effects of cooling on the distribution 
of particles injected at the termination shock are very sim- 



ilar to those calculated by Kardashev (1962) under the 
assumption that each particle spends the same length of 
time within the cooling zone. The upper cut-off moves to- 
wards lower Lorentz factors as r rises, but remains quite 
sharp. A low energy tail at Lorentz factors less than those 
injected is almost absent for r < 1, where the cooling 
time of the lowest energy injected particle is longer than 
the transport time tdiff- For larger values of r, the parti- 
cles escaping from the head region have a sharply peaked 
distribution centred at 7 w 1/r. 

For diffusive transport (a = 1) the spread in escape 
times is appreciable, resulting in a smoother cut-off and 
earlier development of a low energy tail at Lorentz fac- 
tors 7 < 1. Results for this case have been presented 
by |Wilson (1975|) and [Valtaoja (1984|). In sub-diffusion 



(a < 1) these trends are even stronger. Even for the rela- 
tively fast transport through the cooling region implied by 
t = 10 -2 , the particle spectrum is steepened significantly 
over the entire injected range from 7 = 1 to 7 = 10 4 . At 
low Lorentz factors a power-law tail is noticeable already 
at r = 1CP 4 . The very large spread in escape times implied 
by this model can be seen from the distribution plotted for 
t = 100. A substantial number of particles leave the head 
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region with Lorentz factors of 7 = 10. These particles have 
escaped within 10 -3 of the characteristic transport time 
tdiS- On the other hand, the distribution for r = 10 
reveals particles escaping with 7 = 0.2, which have been 
retained in the head region for over 10 4 x t^is- 

In summary, the main effect of cooling in the high-loss 
head region on the particle distribution is not to steepen 
it into a new power-law, but to impose a cut-off of varying 
degrees of sharpness, according to the transport proper- 
ties. The adiabatic expansion losses suffered by these par- 
ticles before entering the lobes shifts the escaping spec- 
trum in energy and amplitude, but does not change its 
form. 

4.2. The lobe 

The effect of synchrotron and inverse Compton losses that 
are constant in time on a power-law distribution Qiobc oc 
7 _s of injected electrons is well known. However, losses 
in combination with time-dependent adiabatic expansion 
and continuous injection — processes present in Eq. ( pi]) 
- considerably complicate this simple picture. Analytic 
solutions can be found for the regimes in which one or 
other of these processes dominates. Two of these are given 
in the Appendix and are briefly discussed in this section. 
In each case it is assumed that a power-law distribution 
is continuously injected in the range 71 < 7 < 72. 

4.2.1. Adiabatic losses 

In this case Eq. ( p^ ) becomes 7 = — 67/t. The solution of 
Eq. (|2|) is given by Eq. (p9| ) . It is characterised by a range 
over which the distribution has the same slope as the in- 
jection function N oc 7~ s and accumulates approximately 
linearly with time. Above a critical value of the Lorentz 
factor, which decreases with time, the distribution cuts off, 
tending to zero at 72. Below 71, a power-law distribution 
is established with A oc —1 + 1/6 which also accumulates 
approximately linearly with time. (In the cases of rele- 
vance here, 6 as 0.4.) The distribution vanishes below a 
lower cut-off, that decreases with time. 

4.2.2. Inverse Compton losses 

In this case Eq. becomes 7 = — &i C 7 2 , and the solution 
of Eq. (|2l]) is given in Eq. (|30|). Below a time-dependent 
break frequency the distribution has the same slope as the 
injection function N cx 7~ s and accumulates steadily with 
time. Above the break frequency the distribution steepens 
by unity: N oc 7 _s_1 and remains constant in time, finally 
cutting off at the maximum injected Lorentz factor 72. For 
7 < 71 a cooling distribution N cx 7~ 2 is established. 

Examples of these properties can be seen in Fig. ^|, 
which shows the evolution of the normalised electron dis- 
tribution, A (7, t) for the case of a high jet-power, high 
redshift radio source. The continuous lines represent the 
numerical results and the crosses (where available) the an- 
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Fig. 3. The evolution of the electron distribution in the 
lobe of a high jet power high redshift radio source, in 
the case of a power-law injection spectrum from the 
hotspot. From bottom to top: t = 1, 10 and 200 Myr, 
respectively. Continuous lines represent the numerical re- 
sults and crosses the analytical solutions as described by 
Eqs. © and @. 

alytical solutions, as described in Eq. (|29|) and (|30]). The 
values of the parameters used in this case are: jet power 
Qo = 1.3 x 10 40 W, source redshift z = 2, hotspot radius 
R = 2.5 kpc, 71 = 10 and 72 = 10 7 . The timescale asso- 
ciated with these parameters is to = 2.5 x 10 5 yr and the 
source activity is assumed to start at i s tart = £0; our nu- 
merical results are insensitive to this choice. The efficiency 
was chosen to be rj = 10 -1 . 

The amplitude of the energy spectrum increases with 
time since the rate at which particles enter the lobe is 
higher than the rate at which they cool to lower energies. 
At epochs earlier than (6 s /6ic) 1/ ' a = i C rit ~ 21 Myr syn- 
chrotron losses dominate over Compton losses at all en- 
ergies. Adiabatic losses dominate below a critical Lorentz 
factor given by 7 crit = bt a ~ l /b s » 10 2 [i(Myr)] 057 at early 
times (t < t crit ) and by 7 crit = b/b ic t w 10 4 /i(Myr) at late 
times (t > t cr it). Thus, in Fig. ||, the distribution at 1 Myr 
is dominated by adiabatic losses for 7 < 10 2 and by (time 
dependent) synchrotron losses for 7 > 10 2 . At epochs 
100 Myr and 200 Myr Compton losses dominate above 
7 = 100 and 50 respectively. Because the synchrotron 
losses decrease with time, even the loss-dominated part of 
the distribution grows significantly between 1 and 10 Myr. 
When Compton losses arc decisive, the loss dominated 
part remains almost constant. Below 7 = 10 (the mini- 
mum Lorentz factor at which injection occurs) a power- 
law distribution of index — 1 + 1/6 ~ 1.6 is established by 
adiabatic losses. 
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5. Results 

In this section we present some results concerning the 
spectral evolution of the lobes of classical double radio 
sources, derived by integrating the electron energy distri- 
bution in the lobes, N(j,t), over the electron emissivity 
given in Eq. (^7j) . The electron energy distribution is com- 
puted numerically by integrating Eq. (|24|) using the dif- 
ferential rate at which electrons are injected into the lobe 
given in Eqs. © and @. 

The parameters we use to specify a radio source are: 

1. The jet power Qq 

2. The efficiency r\ (< 1) with which particles are accel- 
erated at the terminal shock front, which we set to 
unity in case A (without re-acceleration), and to 0.4 
in case B (with re-acceleration). The specific power is 
linearly proportional to r\. 

3. The area A of the terminal hot spot — we take A = 
7T x (2.5) 2 kpc 2 . 

4. The redshift z of the source 

5. The density of the surrounding matter at a fiducial 
distance, according to Eq. (|J). 

6. The density profile, for which we take [3 — 1.5. 

7. Th 



10. 



puwei-law index p of particles accelerated al the 
teri'niiial shuck. Depending uii llie Luienlz factor uf Lhe 

jet, one expects p w 2 to 2.3. We adopt the value 2.23 
( [Kirk et al, 2000| ). 

The minimum Lorentz factor of the accelerated parti- 
cles. This should be somewhat larger than the Lorentz 
factor of the jet. The results are not sensitive to its 
value for the relevant range of spectral indices. 
The ratio t of the typical transport time idiff to cooling 
time i CO oi at 7 = 7 m i n . This is essentially a free param- 
eter. Although spherical geometry is used to motivate 
the transport model in Sect. |^, the resulting concept 
of a stochastic spread in particle escape times from the 
high magnetic field head region is much more generic. 
However, as in the case of standard diffusion, the model 
assumes infinite velocity of a fluid element within the 
head region. This could lead to unphysical results if 
idiff is chosen to be too small. In our case, if i? eS c is the 
minimum distance a fluid element must travel before 
exiting the head region then the ratio of escape time 
to t coo \ must exceed 10~ 4 7 m i n ((5/10 40 W)(i? esc /1 kpc) 
for all particles. In the case of sub-diffusion, where the 
individual escape times have a large spread about the 
mean, r must exceed this value by a factor of a few. 
Thus, in our computations we take r > 2 x 10 -3 . 
The maximum Lorentz factor 7 max to which particles 



are accelerated at the terminal shock. In accordance 



with X-ray observations of hot spots we assume 7 max 
10 7 . As a result, r ^ 7min/7max and the results are 
insensitive to this parameter. 
11. The index a that characterises the transport process 
through the head region: < a < 2. We consider three 
values: sub-diffusion a = 0.5, diffusion a = 1 and supra 
diffusion a = 1.5. 
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Fig. 4. Tracks for two high redshift, high power sources 
for case A (no re-acceleration) , with paramet ers as given in 
Table The 'standard' track corresponds to Blundcll ct a] 
(1999|) , Fig. 13. On each track the dots indicate the source 



age, starting at 1 Myr and uniformly spaced in intervals of 
5Myr. The solid line shows the track for diffusion (a = 1) 
the dashed line for sub-diffusion (a = 0.5) and the dotted 
line for supra-diffusion (a = 1.5). Crosses and open circles 
indicate FR II sources of the 3CR revised sample, with 
z > 0.5 and z < 0.5, respectively. 

In addition we distinguish two cases: 

— Case A: particles suffer time-dependent adiabatic and 
radiative losses on passing from the termination shock 
to the lobes. 

— Case B: adiabatic losses between the termination shock 
and the lobes are compensated by re-acceleration in 
the head region. In this case k(t) = 1 in Eq. (20). 



We consider first tracks on the power versus linear size 
diagram of high power, high redshift sources. Figure |i| 
shows two such tracks for case A — no re-acceleration 
of particles after the termination shock — for parame- 
ters given in Table [|. Three different transport models 
are shown (a = 0.5, 1, and 1.5) and we have chosen 
r = 2 x 10~ 3 . This figure also shows the positions of FR II 
sources in the 3CR revised catalogue, (Laing et al, 1983) 



according to the classifications by Jackson & Rawlinge 



(1997D , as updated a nd used by |Willott et al (1999| ) and 
published on-line by |Willott (2002| ). The linear size was 
computed assuming the cosmological parameters qo = 0, 
Hq = 50kms _1 Mpc -1 . Sources with redshift z > 0.5 are 
indicated by a cross, those with z < 0.5 by an open circle. 

The first track, lab elled 'standard' can b e directly com- 
pared with Fig. 13 of Blundcll ct al (1999 ), to reveal the 
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new features of our model. We find a much sharper decline 
in the source power with increasing radius than predicted 
by Blundcll ct al (1999| ) . For these computations we adopt 
the minimum plausible value of the transport time: r = 
2 x 1CP 3 — increasing r always leads to an even steeper de- 
cline. The track passes well below the large {D > 600 kpc) 
sources with specific power P151 > 3 x 10 28 WHz _1 and 
decays to P 15 i « 7x 10 26 W Hz" 1 within less than 26 Myr. 
The reason for this behaviour is the lack of particles of high 
Lorentz factor entering the lobes. Synchrotron losses in the 
head region effect a cut-off of the injected distribution in 
our model of acceleration and transport. After suffering 
the strong expansion losses associated with the pressure 
difference between primary hot spot and lobe, this cut- 
off affects the number of particles available to radiate at 
151 MHz, which, as a result, declines steeply with age. In 
order to account for the larger sources, it is necessary not 
only to assume a much larger power Qo as 6.5 x 10 40 W, but 
also a much higher surrounding density p ~ 10 _22 kgm~ 3 
at 10 kpc. Such a track is shown in Fig. [| labelled 'high 
density'. At age roughly 22 Myr this track coincides with 
the largest high power source in the 3CR revised cata- 
logue. If sources existed with these extreme parameters 
they should also be visible at much higher specific pow- 
ers, when less than 10 Myr old. Such objects, however, do 
not appear in the catalogue. We therefore conclude that 
case A tracks, in which particle acceleration occurs only 
at the termination shock are in conflict with the data. 

This conclusion is unaffected by the type of transport 
present in the head region. The dashed lines in Fig. [| show 
the effects of sub-diffusion and the dotted lines that of 
supra-diffusion on the 'standard' and 'high density' tracks. 
For ages greater than 5 Myr, sub-diffusion always results 
in greater specific power than diffusion (a = 1), whereas 
supra-diffusion gives a lower specific power. This arises 
because sub-diffusion produces a greater spread in escape 
times from the head region than does diffusion. The upper 
cut-off in the Lorentz factor of particles entering the lobes 
is correspondingly higher, whereas the opposite is the case 
for supra-diffusion. Nevertheless, the sub-diffusion track 
for the standard source is unable to account for the largest 
high redshift sources, and the supra-diffusion track for a 
high density source would predict more medium-sized high 
specific power sources than are observed. 

Turning to case B, in which the adiabatic losses be- 
tween head region and lobe are compensated by a re- 
acceleration process operating downstream of the termi- 
nation shock, we present in Fig. [5] two tracks in the P—D 
diagram, one corresponding to a high redshift, high power 
source and one to a low redshift, low power source. Once 
again, the parameters are chosen to enable direct compar- 
ison with the work of ( Blundcll et al, 199S| ) and are shown 
in Table [l]. The value of r is the same as in Fig. [|. 

A prominent feature of these tracks is the steepening 
as the source ages. To a lesser extent, this is also seen in 
case A, Fig. an d is p resent in the tracks computed by 
Kaiser et al (1997 ) and Blundcll ct al (199S ). It arises at 
the point where inverse Compton losses become impor- 



Table 1. Parameters of the models displayed in Figs, 
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Fig. 5. Tracks for a high redshift, high power source 
and a low redshift, low power source for case B (includ- 
ing re-acceleration), with parameters as given in Table 0. 
Diffusive transport a — 1 is shown by a solid line, on which 
the source age is indicated by the dots, starting, for the 
high z source, at 1 Myr and, for the low-z source, at 4 Myr 
and uniformly spaced in intervals of 5 Myr in each case. 
The dashed lines indicate sub-diffusion a = 0.5 and the 
dotted lines supra-diffusion (a — 1.5) Crosses and open 
circles indicate FR II sources of the 3CR revised sample, 
with z > 0.5 and z < 0.5, respectively. 

tant in the lobes. Whereas the synchrotron losses continu- 
ally weaken with time, the inverse Compton losses remain 
constant, causing a more dramatic fall off in the specific 
power. 

Because of the dependence of the source size on its 
power (see Eq. (^)) the low z source in Fig. || is sys- 
tematically older at a fixed source size than is the high 
z source — a potentially important fact in the interpre- 
tation of the correlation between source size and spec- 
tral index. In Fig. ^ we plot the particle distribution in 
the lobes, multiplied by 7 P , against Lorentz factor 7 for 
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Fig. 6. The particle distribution function in the lobes, 
multiplied by 7 P , where p is the index of the injected dis- 
tribution for the high and low redshift sources shown in 
Fig. [jj. The solid lines depict the distribution for the high 
z source, the dashed lines for the low z source. For each 
source the distribution is plotted when it achieves the sizes 
56, 213 and 1127 kpc. For fixed size, the higher redshift 
source has the softer spectrum. For fixed redshift, the older 
(larger) source is softer. 



both the high z source (solid line) and the low z source 
(dashed line). Despite the age difference, it can be seen 
that where cooling steepens the spectrum, the high power, 
high z source always has a spectrum that is softer than 
that of the low power, low z source. Furthermore, for fixed 
redshift, the source spectrum systematically softens with 
increasing size, and, therefore, also with increasing age. 



of escape times for the particles from the high-loss head 
region. Compared to the case of standard diffusion, these 
times have a much wider distribution in sub-diffusion and 
a much tighter one in supra-diffusion. The new formalism 
also requires a parameter r that corresponds to the av- 
erage residence time in units of the cooling time. Thus, 
in total, one more parameter is used than the standard 
Kardashcv (1962] ) model with a single escape time for 
all particles and one parameter fewer i s needed than the 



model introduced by [Blundefl et al (1999|) , which requires 
two break frequencies and a 'leakiness'. 

The self-similar description of the hydrodynamical evo- 
lution of classical double radio sources implies that the 
lobe pressure falls with time. Because the pressure in the 
primary hot-spots does not appear to undergo this evolu- 
tion, the particles injected at there in older sources have to 
overcome ever more severe adiabatic losses. These losses 



were not modelled by Kaiser et al (1997 ), wh o specified the 
distrib ution of particles entering the lobes. Blundcll et a] 
(1999| ) pointed out their importance and included them in 



their P-D tracks. However, because the transport model 
did not include an upper cut-off of injected particles, but 
merely a spectral break, they were still able to account 
for the observed population of sources. We find that the 
adiabatic losses between primary hot spot and lobe pre- 
vent this when our transport model is adopted. The large, 
high luminosity sources can be explained only as very pow- 
erful jets in very dense environments, in which case they 
should have smaller, brighter predecessors, that are absent 
from the 3CR revised catalogue. There could be several 
reasons for this. One possibility is that at high redshift, 
small sources with powerful jets exist, but are not bright 
in the radio because of absorption. The missing progeni- 
tors could then be powerful, high red shift analogues of th e 
"Gigahertz Peak Spectrum" sources ( Bicknell et al, 1997 ). 
Another is that the jet power Qq remains quite weak dur- 
ing the early evolution of a high-redshift source, becom- 
ing strong only after ~ 10 Myr. Intermittent jet activity 
in FR II sources has already been proposed ( Reynolds & 



6. Conclusions 

We have introduced a new description of the accelera- 
tion and transport of particles in the hot spots and lobes 
of FR II radio galaxies. Using the standard approach of 
assuming a power-law spectrum injected at the termina- 
tion shock, we model the subsequent transport through 
the head region into the lobes using a formalism which 
permits the investigation of 'anomalous' regimes. These 
transport regimes are described by the single parameter 



a, whi:h determines the time-dependence of the mean- 
square displacement of a particle: (Ar 2 ) cx t a . Standard 
diffusion, in which the flux is proportional to the gra- 
dient of the particle density (Fick's law) corresponds to 
a = 1 . The anomalous regimes of sub-diffusion and supra- 
diffusion correspond respectively to a < 1 and a > 1. In 
the current application, the most important physical dif- 
ference between these regimes is that they permit a range 



Begelman, 1997 ), but as a solution to the opposite problem 
of an observed cwer-abundance of small powerful sources 
(at lower redshift). In general, it seems easier to imag- 
ine the jet power to be larger during earlier phases of a 
source's life, leading to a faster decay of the specific source 
power with time. 

However, the solution we propose to this problem is 
that acceleration is not confined to the primary hot spot, 
but occurs throughout the head region. Independent evi- 
dence from individual sources in the infra-red, optical and 
X-ray bands (|Mciscnhcimcr et al, 1996 , 1997; Perley et al 



1997; Wagner & Krawczynski, 2000) supports this idea. A 



plausible mechanism is repeated encounters with weaker 
shock fronts which may permeate this turbulent region. 
We have investigated the P-D tracks and spectra of a 
model in which it is assumed that this distributed acceler- 
ation compensates the expansion losses between primary 
hot spot and lobes, without significantly modifying the 
spectrum. Our tracks are more akin to those presented by 
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Kaiser et al (1997) and wc have investigated their sensi- 
tivity to the new transport parameter a. We have shown 
that some of the properties noted by Blundcll ct al (1999| ) , 



such as a redshift /spectral index and a size/spectra index 
correlation also emerge from our model. 
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7. Appendix 



Analytic solutions of Eq. (J21j) are presented for the cases in which either only adiabatic losses or only synchrotron 
losses are present, assuming an injection function given by 



/lobe 



(7) = 



(?o7 s for 71 < 7 < 72 



otherwise 



7.1. Adiabatic losses 

The solution of Eq. [H], provided that 71 < ^{t / tstiat)~ b , is: 

r 0, 



<?0* S tart7 



-(*)*(.£) *-[(*) (rib) 



/( 7 , t) = ^ 90*^7 ° ( _^_) 

A y Estart / \ *start / 



7 < 7ai 



, 7 ai < 7 < 7i 



7i < 7 < 7q 2 



O0t S tart7 



0, 



\ *start / 



1- ( J - 

72 



7 > 72 



-6 



with A = b(s - 1) + 1 and 7 0l = 7i(t/i start ) b , 7 Q2 = 72(*/*start) 
The solution of Eq. |l], for 71 < [7.7 1 + 6 ic (t - i st art)]~\ is: 

f 0, 7 < 7 Cl 

^7- (s+1) {(^) s_1 [i-&ic7(*-wt: 

/(7, *) = { ^ 7 -( s +D{l - [1 - b ic7 (t - istart)]^ 1 } 
J^ 7 -(s + l) 1 _ ,_ V - I 



, 7cj < 7 < 7i 



1- (-) 

v 72 / 



0. 



7i < 7 < 7c 2 
7c 2 < 7 < 72 
7 > 72 



(28) 



(29) 



(30) 



where K = q /b ic (s - 1), 7 C1 = [7i 1 + ho(t - t staxt )] 1 and j C2 = [j 2 1 + 6 ic (t - i start )] . 

In the case that adiabatic losses determine the spectrum of a source the solution of Eq. (p|) for 71 > j a2 , is: 



r 0, 



OQ*Btart7 
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1 - 



71 / V *et 



/( 7 ,t) = | ^f-»» (^) 7l -^[l-(^) 
?ot s t art7 ~ p f t A i_(sl) 

A ytstart / \72 J 



A/6' 



7 < 7a! 



7ai < 7 < 7a 2 



7a 2 < 7 < 71 



7i < 7 < 72 



(31) 



0, 



7 > 72 
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In the case that inverse compton losses determine the spectrum of a source the solution of Eq. (|l|) for 7 > 7 C2 , is: 

r 0, 7 < 7ci 

^ 7 - 2 { 7l 1 - p -[ 7 - 1 -6 ic (t-t s 1 





( 7 \ v ~ x ~ 


1 - 







}, 7ci < 7 < 7c 2 
7c 2 < 7 < 7i 

7i < 7 < 72 

7 > 72 



(32) 
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